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New  thermal  rate  equations  were  developed  by  taking  the  temperature  dependences  of  the  electrical 
resistivity  p  and  thermal  conductivity  k  of  the  thermoelectric  (TE)  materials  into  the  thermal  rate  equa¬ 
tions  on  the  assumption  that  they  vary  linearly  with  temperature  T.  The  relative  energy  conversion  effi¬ 
ciency  rjlrjo  for  a  single  TE  element  was  formulated  by  approximate  analysis,  where  rj  and  rj0  are  the 
energy  conversion  efficiencies  derived  from  the  new  and  conventional  thermal  rate  equations,  respec¬ 
tively.  Applying  it  to  Si-Ge  alloys,  the  temperature  dependence  of  p  is  stronger  than  that  of  k,  so  the  for¬ 
mer  has  a  more  significant  effect  on  rj/rjo  than  the  latter.  However,  the  degree  of  contribution  from  both  of 
them  to  rjho  was  a  little  lower  than  1%  at  the  temperature  difference  AT  of  600  K.  When  the  temperature 
dependence  of  k  was  increased  to  become  equal  to  that  of  p ,  however,  it  was  found  that  rj/rjo  is  increased 
by  about  10%  at  AT  =  600  K.  It  is  clarified  here  that  the  temperature  dependences  of  p  and  k  are  also 
important  factors  for  an  improvement  in  rj. 

©  2008  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

There  has  been  considerable  interest  during  the  past  ten  years 
in  finding  new  materials  and  structures  for  use  in  clear,  highly  effi¬ 
cient  cooling  and  energy  conversion  systems  [1,2].  The  thermo¬ 
electric  (TE)  figure  of  merit  ZT  provides  a  measure  of  the  quality 
of  such  materials  for  applications  and  is  defined  by  ZT  =  a2T/p?c, 
where  a  is  the  Seebeck  coefficient,  p  is  the  electrical  resistivity,  k 
is  the  thermal  conductivity  and  T  is  the  absolute  temperature  [1]. 
The  increase  in  ZT  leads  directly  to  improvement  in  the  energy 
conversion  efficiency  of  TE  generators  and  in  the  cooling  efficiency 
of  Peltier  modules  [3].  Much  effort  has  been  made  to  raise  the  ZT  of 
TE  bulk  materials  for  energy  conversion  efficiency,  so  there  have 
been  some  improvements  in  ZT  (for  example,  ~2  at  700  K  for  the 
p-type  TAGS  (GeTe-AgSbTe2)  [4],  about  1.4  at  620  K  for  the  p-type 
CeFe3.5Coo.5Sbi2  skutterudite  compound  [5]  and  about  1.5  at  700  K 
for  the  p-type  (Zr-Hf)NiSn  half-Heusler  compound  [6]).  However, 
even  their  values  are  not  sufficient  to  improve  dramatically  the  en¬ 
ergy  conversion  efficiency  ( rj ). 

The  energy  conversion  efficiency  rj  has  been  calculated  analyti¬ 
cally  using  the  well  known  thermal  rate  equations  at  the  cold  and 
hot  junctions  of  a  TE  element  [3].  The  change  in  a  with  tempera¬ 
ture  has  been  taken  into  the  conventional  thermal  rate  equations. 
However,  the  variations  in  p  and  k  with  temperature  were  ne¬ 
glected,  although  generally  they  have  different  values  at  the  hot 
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and  cold  junctions  [7].  The  traditional  expressions  for  t]  obtained 
from  the  conventional  thermal  rate  equations  should,  thus,  not 
be  exactly  correct.  For  this  reason,  rj  has  been  calculated  numeri¬ 
cally  by  solving  the  differential  equations  governing  TE  properties 
by  the  software  tool  in  which  the  correct  temperature  dependence 
of  TE  properties  is  taken  into  account  [8].  However,  such  a  calcula¬ 
tion  has  never  clarified  analytically  the  degree  to  which  the  tem¬ 
perature  dependence  of  each  TE  property  has  an  effect  on  rj.  To 
make  it  clear,  we  made  an  attempt  to  take  the  difference  between 
the  TE  properties  at  the  hot  and  cold  sides  of  a  TE  element  into  the 
thermal  rate  equations  on  the  assumption  that  all  of  the  TE  prop¬ 
erties  change  linearly  with  temperature.  The  new  thermal  rate 
equations  were,  thus,  formulated  by  taking  the  new  temperature 
dependences  of  both  p  and  k  into  the  thermal  rate  equations. 
The  reason  that  only  the  linear  temperature  dependence  of  TE 
properties  was  employed  here  is  simply  to  make  clear  the  effect 
of  the  temperature  dependence  of  TE  properties  on  rj.  It  is  clarified 
for  Si-Ge  alloys  [9]  in  this  paper  how  rjlrj0  varies  with  changes  in 
the  sign  and  magnitude  in  the  gradients  of  p  and  k  for  T,  where 
rj/rjo  is  the  relative  energy  conversion  efficiency  calculated  using 
the  new  and  conventional  thermal  rate  equations,  respectively. 
The  reason  that  Si-Ge  alloys  were  employed  for  the  calculation 
of  f]lrj0  is  that  their  TE  properties  indeed  change  linearly  with  tem¬ 
perature  in  the  broad  temperature  region  from  400  to  1000K.  rj 
was  then  derived  analytically  on  the  assumption  that  the  terms 
dependent  on  the  temperature  difference  AT  imposed  on  a  TE  ele¬ 
ment  are  sufficiently  small  compared  with  those  independent  of 
AT. 
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Nomenclature 

TE 

thermoelectric 

D 

normalized  temperature  dependence  of  Z(T)  at  T  (de- 

Z 

thermoelectric  figure  of  merit  (K_1) 

fined  as  D  =  (dZ(T)/dT))/Z(T))  (K_1) 

ZT 

dimensionless  thermoelectric  figure  of  merit 

Qh 

heat  energy  at  hot  junction  of  TE  element  (W) 

T 

mean  temperature  of  TE  element  (given  by  T=  (Th+Tc/2) 

Qc 

heat  energy  at  cold  junction  of  TE  element  (W) 

(K) 

Q 

difference  (Qh  -  Qc)  in  heat  energy  between  both  end 

AT 

temperature  difference  between  both  ends  of  TE 

junctions  of  TE  element  (W) 

element  (K) 

S 

cross  sectional  area  of  TE  element  (m2) 

Tz 

temperature  at  position  z  of  TE  element  (K) 

L 

length  of  thermoelectric  material  of  TE  element  (m) 

Th 

temperature  at  hot  junction  of  TE  element  (given  by  T 

Rl 

load  resistance  (Q) 

+AT/2)  (K) 

I 

electric  current  through  TE  element  (A) 

Tc 

temperature  at  cold  junction  of  TE  element  (given  by 

T— AT/2)  (K) 

Greek  letters 

m 

thermal  conductance  given  by  I<(T)  =  Sk(T)IL  (W/K) 

a 

Seebeck  coefficient  of  material  (V/K) 

R(T) 

electric  resistance  given  by  R(T)  =  Lp(T)/5  (Q) 

a(Th) 

Seebeck  coefficient  at  hot  junction  of  TE  element  (V/K) 

<K(Th)) 

thermal  conductance  averaged  over  half  (hot)  side  of  TE 

oc(Tc) 

Seebeck  coefficient  at  cold  junction  of  TE  element  (V/K) 

element  (W/I<) 

K 

thermal  conductivity  of  material  (W/mK) 

(I<(Tc)) 

thermal  conductance  averaged  over  half  (cold)  side  of 

P 

electrical  resistivity  of  material  (Q  m) 

TE  element  (W/I<) 

y 

parameter  defined  as  (1  +ZT)- 1/2 

(R(Th)) 

electrical  resistivity  averaged  over  half  (hot)  side  of  TE 

P(Th) 

Thomson  coefficient  at  hot  junction  of  TE  element 

element  (Q) 

(defined  as  Th(da/dT))  (V/K) 

ro> 

electrical  resistivity  averaged  over  half  (cold)  side  of  TE 

P(Jc) 

Thomson  coefficient  at  cold  junction  of  TE  element 

element  (Q) 

(defined  as  Tc(da/dT))  (V/K) 

A 

normalized  temperature  dependence  of  a(T)  at  T 

n 

energy  conversion  efficiency  for  B  /  0  K-1  and  CV  0  K_1 

(defined  as  A  =  (da(T)/dT))/a(T))  (IT1) 

(%) 

B 

normalized  temperature  dependence  of  p{T)  at  T 
(defined  as  B  =  (dp(T)/dT))/p(T))  (IT1) 

ho 

energy  conversion  efficiency  for  B  =  C  =  0  K  1  (%) 

C 

normalized  temperature  dependence  of  k(T)  at  T 
(defined  as  C=  (dK(T)/dT))/K(T))  (IT1) 

The  purpose  of  this  paper  is  to  investigate  to  what  extent  the 
temperature  dependences  of  p  and  k  of  TE  materials  have  an  effect 
on  rj  of  a  single  TE  element  and  to  show  whether  rj/rj0  of  a  TE  ele¬ 
ment  can  be  enhanced  significantly  by  improving  the  temperature 
dependences  of  p  and  k. 

2.  Analysis 

2.1.  Temperature  dependence  of  thermoelectric  properties 

Let  us  consider  a  single  TE  element  in  which  a  TE  parallelepiped 
is  welded  with  two  electrically  conducting  materials,  as  shown  in 
Fig.  1.  When  a  temperature  difference  is  imposed  on  a  TE  element, 
a  voltage  is  generated  between  both  ends  of  a  TE  element,  and  a 
current  flows  through  a  circuit.  The  current  can  then  be  used  to 


/ 


Fig.  1.  A  single  TE  element  (possessing  the  length  L  and  cross-sectional  area  S) 
sandwiched  between  two  metals  on  which  the  temperature  difference  AT  is 
imposed,  where  Th  and  Tc  are  the  temperatures  at  the  hot  and  cold  junctions,  T  is  the 
mean  element  temperature  and  /  is  a  current  flowing  through  a  circuit  composed  of 
a  TE  element  and  a  load  resistance  RL. 


power  a  load,  thus  converting  the  thermal  energy  into  electric  en¬ 
ergy.  For  a  majority  of  cases,  a  one  dimensional  model  has  been 
sufficient  to  obtain  quantitatively  correct  performance  data  [10]. 
Assuming  that  the  temperature  difference  between  the  hot  (Th) 
and  cold  (Tc)  junctions  of  a  TE  element  is  AT  and  the  mean  temper¬ 
ature  is  T,  Th  and  Tc  are  expressed  as  T  +  AT/2  and  T  -  AT/2,  respec¬ 
tively.  Thus,  the  energy  conversion  efficiency  rj  should  depend 
strongly  on  the  temperature  dependence  of  the  TE  properties  of 
TE  materials  because  the  TE  properties  are  different  at  the  hot 
and  cold  junctions. 

Here,  let  us  assume  that  the  Seebeck  coefficient  a(T),  the  electric 
resistance  R(T )  and  the  thermal  conductance  K(T)  of  a  TE  element 
change  linearly  with  an  increase  of  T  over  a  working  temperature 
range,  where  T  is  the  mean  temperature.  a(Th)  at  the  hot  junction 
and  a (Tc)  at  the  cold  junction  of  a  TE  element  are  then  expressed 
as  a  function  of  AT  as 

a(T±AT/2)  =  a(T)(l  ±AAT/2),  (1) 

and  R(Th )  and  I<(Th )  at  the  hot  side  and  R(TC )  and  K(TC)  at  the  cold 
side  of  a  TE  element  are  given  in  the  same  way  by 

R(T  ±  AT/2)  =  R(T)(  1  ±  BAT/2)  (2) 

and 

K(T  ±  AT/2)  =  /C(T)(1  ±  CAT/2),  (3) 

where  A  =  (da(T)/dT)/a(T),  B  =  (dR(T)/dT)/R(T)  and  C=  (d/C(T)/dT)/ 
K(T).  The  difference  in  a  between  both  junctions  has  been  taken  into 
the  conventional  thermal  rate  equations,  but  the  differences  in  R 
and  I<  between  the  hot  and  cold  sides  have  always  been  neglected. 
In  the  conventional  thermal  rate  equation,  a  has  been  replaced  with 
an  average  value  (a)  over  a  working  temperature  range,  and  the 
resistances  of  the  hot  and  cold  sides  were  replaced  by  R/2  so  that 
the  Joule  heat  is  divided  equally  between  the  hot  and  cold  sides. 
As  for  K,  the  K  value  at  the  mean  temperature  has  hitherto  been  em- 
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ployed  for  the  thermal  rate  equations,  so  that  I<(T)  was  regarded  as  a 
constant  throughout  a  TE  element.  Thus,  the  temperature  depen¬ 
dence  of  R(T )  and  I<(T )  has  never  been  taken  into  account  in  the  con¬ 
ventional  thermal  rate  equations. 

When  p(T)  and  k(T )  change  linearly  with  T,  it  is  shown  that  if  we 
divide  R(T )  of  a  TE  element  into  ( R(Th ))  and  ( R(TC ))  averaged  over  a 
half  (hot  or  cold)  side  of  a  TE  element  and  replace  I<{T)  with  (I<(T)) 
averaged  all  over  a  TE  element,  they  can  be  expressed  using  the 
average  gradients  in  p  and  k  for  T.  The  reason  that  the  resultant 
(K(T))  was  employed  here  for  the  thermal  rate  equations  instead 
of  (I<{Th))  and  (I<{TC))  is  that  the  thermal  conductance  must  be  con¬ 
sidered  as  a  total  thermal  conductance.  As  described  in  the  Appen¬ 
dix,  when  the  temperature  within  a  TE  parallelepiped  changes 
linearly  along  its  length  direction  (see  Fig.  1),  the  average  (B(Th)) 
and  (B(TC))  of  the  hot  and  cold  sides  are  then  expressed  as 

(R(Th))  =  (K(T)/2)(1+ BAT/4)  (4) 


By  substituting  Eqs.  (1),  (4),  (5),  (9)  and  (10)  into  Eqs.  (7)  and 
(8),  the  new  thermal  rate  equations  are  expressed  as 

Qh  =  a(T)ThI-R(T)(  1  +BAT/4)I2/2  +  /C(T)[1  -  (CAT)2/12]AT 

(11) 

and 

Qc  =  oc(T)TcI  +  R(T)(1  -  BAT/4)I2/2  +  K(T)[  1  -  (CAT)2/U]AT, 

(12) 

at  least  when  there  is  no  heat  transfer  to  or  from  the  surroundings 
other  than  a  TE  element.  After  all,  the  term  A  is  canceled  by  that  of 
the  Thomson  coefficient,  as  already  pointed  out  by  Ioffe  [7],  so  that 
only  the  terms  B  and  C2  remain  in  the  thermal  rate  equations.  If 
B  =  C  =  0  K_1,  of  course,  these  thermal  rate  equations  reduce  to  the 
conventional  thermal  rate  equations  [7].  The  energy  difference  Q 
between  Qh  and  Qc  is  given  by 


and 

(R(Tc))  =  (R(T)/2)(l-BAT/4)  (5) 

as  a  function  of  AT.  However,  the  sum  of  (R(Th))  and  (B(TC))  is  equal 
to  R(T).  The  conventional  theory  has  taken  B  =  OK-1,  but  we  will 
take  here  B  ±  0  Kr1  because  BAT/4  has  a  value  of  0.19  on  average 
at  AT  =  600  I<  in  Si-Ge  alloys  [8],  a  value  that  is  not  negligibly  small 
compared  with  1 . 

When  AT  is  applied  to  a  TE  element,  the  resultant  ( I<(T ))  aver¬ 
aged  all  over  a  TE  element  is  expressed  as 


Q  =  Qh  -  Qc  =  0C(T)AT  I  -  R(T)I2  =  RlI2,  (13) 

where  RLI2  is  the  electric  energy  consumed  in  a  load  resistance  RL 
[7].  From  Eq.  (13),  the  current  /  flowing  through  a  circuit  is  obtained 
simply  as 


a(T)AT 

W+ WT) 


(14) 


The  energy  conversion  efficiency  p  of  a  TE  element  is  defined  as  QJ 
Qh  [7]  and  is  expressed  as 


(K(T))=K(T)[1-(CAT)2/12],  (6) 

as  described  in  the  Appendix,  if  (CAT)2/12  is  very  small  compared 
with  1,  where  I<(T)  =  Sk{T)IL.  The  resultant  (K(T))  depends  on  the 
term  of  C2,  not  on  C.  It  is,  thus,  found  that  the  resultant  (K(T))  tends 
to  decrease  significantly  with  an  increase  of  AT,  irrespective  of  the 
sign  of  C. 


n  =  Q/Qh 

= _ Rif _ 

a(T)Th/-K(T)(l  +BAT/4)/2/2  +  /C(T)[l  -  (CAT)2/12]AT 

(15) 

By  substituting  Eq.  (14)  into  Eq.  (15)  and  setting  BL/B(T)  =  m,  p  is 
rewritten  as 


2.2.  p  derived  by  taking  into  account  the  temperature  dependence  of 
thermoelectric  properties 

Here,  we  make  an  attempt  to  take  the  temperature  dependenc¬ 
es  of  TE  properties  into  the  energy  conversion  efficiency  p.  First  of 
all,  we  derive  the  new  thermal  rate  equations  using  a(Th),  a(Tc), 
(B(Th)),  (R(TC)),  p{Th)  and  /?(TC),  where  p(Th)  and  p(Tc)  are  the  Thom¬ 
son  coefficients  given  by  the  relations  p(Th)  =  Th(da/dT)  and 
p{Tc)  =  Tc(da/dT)  [11].  When  AT  is  applied  to  a  TE  element  and 
the  Seebeck  and  Thomson  effects  result  in  power  generation,  the 
thermal  rate  equations  for  a  TE  element  can  be  expressed  as 

Qh  =  [a(Th)  -  /J(Th)AT/2Th]Th/  -  (R(Th))I2  +  (I((T))AT  (7) 

and 

Qc  =  [a(Tc)  +  P(JC)AT/2TC]TCI  +  ( R(TC))1 2  +  (I((T))AT,  (8) 

where  T  is  the  mean  TE  element  temperature  of  (Th  +  Tc)/2,  (. K(T ))  is 
the  resultant  average  thermal  conductance  under  AT,  and  Qh  and  Qc 
are  the  thermal  energies  heated  and  cooled  at  the  hot  and  cold 
junctions,  respectively  [12].  The  first,  second  and  third  terms  in 
both  equations  represent  the  Peltier  and  Thomson  effects, 
Joule  heating  and  thermal  conduction,  respectively.  When  (B(Th)>  = 
(R(TC))  =  R(T)1 2  and  (K(T))  =  I<{T),  these  equations  correspond  to  the 
thermal  equations  reported  by  Min  et  al.  [12].  Using  the  relation 
between  the  Thomson  coefficient  p  and  A  defined  in  the  previous 
subsection,  the  following  relations  are  obtained  as 

P(Th)  =  ^a(T)  (9) 


n  =  Q/Qh 

ZATm 

~  (1  +  m)2[l  -  (CA7-)2/12]  +ZTh(l  +  m)  -  (1  4-  BAT/4)(ZAr)/2  ’ 

(16) 

where  Z=  a(T)2/B(T)/C(T).  By  making  drjldm  equal  to  zero,  the  qua¬ 
dratic  equation  of  m  for  maximum  r\  is  obtained  as 

(1  -m2)[l  -(CAT)2/12]+ZTh-(l  -  BAT/4) (ZAT)/2  =  0.  (17) 

When  the  values  of  (CAT)2/12  and  BAT/4  are  very  small  compared 
with  1,  Eq.  (17)  is  rewritten  as 

m2  =  y~2[  1  +  (1  -  y2)(AT)2(2TC2  -  3B)/24T]  (18) 

and 

m  =  y-1[l  +  (1  -  y2)(AT)2(2TC2  -  3B)/48T],  (19) 

where 

7  =  (1  +ZT)_1/2.  (20) 

By  substituting  Eqs.  (18)  and  (19)  into  Eq.  (16),  moreover,  the  en¬ 
ergy  conversion  efficiency  p  is  expressed  by  a  similar  approxima¬ 
tion  as 


V  =  ho 


1  + 


4(yTc  +  Th) 


B(i-r)  i  £1(1+ 7) 


(21) 


where  p 0  is  the  original  energy  conversion  efficiency  derived  from 
the  conventional  thermal  rate  equations  and  is  expressed  as  [7] 


and 

P(TC)  =  TcAa(T). 


(l-y)A  T 

/0  (yTc  +  Th) 


(10) 


(22) 
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Although  Eq.  (21)  is  derived  by  an  approximate  analysis,  it  is 
sufficient  to  investigate  the  effect  of  B  and  C  on  rj.  It  is,  thus,  found 
that  when  a,  R  and  I<  have  a  linear  dependence  on  temperature, 
only  the  terms  related  to  the  temperature  dependences  of  R  and 
I<  remain  in  p.  Of  course,  p  is  reduced  to  p0  for  B  =  C  =  0  K-1.  Since 
the  coefficients  of  B  and  C2  are  always  positive  in  sign  because  of 
7(1,  the  increase  in  B  and  C2,  thus,  results  directly  in  an  increase 
in  rj/rjo.  However,  the  ideal  Carnot  efficiency  of  p  is  obtained  in 
the  limit  as  ZT  -►  oo  (7  -►  0),  as  in  the  case  of  p0,  irrespective  of 
the  presence  of  B  and  C. 

The  thermoelectric  figures  of  merit  Z  at  Th  and  Tc  were 
calculated  using  Eqs.  (l)-(3).  When  the  products  of  A,  B  and  C 
and  AT  are  very  small  compared  with  1,  the  expressions  for  Z(Th) 
and  Z(TC)  at  the  hot  and  cold  junctions  are  expressed  approxi¬ 
mately  as 

Z(Th)  =  [a(Th)2/R(Th)K(Th)]  =  Z(T)[  1  +(2A-B-  C)AT/2] 

=  Z(T)(1  +  DAT/2)  (23) 

and 

Z(TC)  =  [a(Tc)2/R(Tc)if(rc)]  =  Z(T)[  1  -  (2 A-B-  QAT/2] 

=  Z(T)(\  -  DAT/2),  (24) 

using  A,  B  and  C  and  AT,  where  Z(T)  =  a(T)2/R(T)J<(T)  and 
D  =  2 A-B-C.  Z(Th)  and  Z(TC)  are,  thus,  found  to  change  significantly 
with  AT  according  to  the  value  of  D.  When  D  =  0  K“\  of  course,  Z 
never  varies  with  changes  in  T.  The  gradient  of  Z  for  T  is  obtained 
from  the  relation  dZ/dT  =  [Z(Th)-Z(Tc)]/AT  as 

dZ/dT  =  ZD  (25) 

using  Eqs.  (23)  and  (24). 

3.  Results  and  discussions 

3A.  Parameters  of  A,  B,  C,  D  and  p  obtained  for  the  p-  and  n-type  Si-Ge 
alloys 

Figs.  2  and  3  show  the  temperature  dependences  of  a,  p ,  k  and  Z 
measured  for  the  p-type  Sio.95Geo.05  alloy  doped  with  B  and  n-type 
Si0.95Ge0.06  alloy  doped  with  P,  which  have  values  of  ZT  =  0.57  and 
0.90  at  1073  K,  respectively  [9].  The  TE  properties  of  Si-Ge  alloys 
change  almost  linearly  with  temperature,  as  shown  in  Figs.  2  and 
3.  For  this  reason,  the  parameters  of  A,  B,  C  and  D  were  calculated 
using  Eqs.  (1),(2),(3)  and  (25)  from  the  experimental  values  of  a,  p, 
k  and  Z  for  Si-Ge  alloys  in  the  temperature  range  from  400  to 
1000  K,  as  listed  in  Table  1.  There  was  no  significant  difference  in 
A,  B ,  C  and  D  between  the  p-  and  n-type  Si-Ge  alloys.  In  both  the 
p-  and  n-type  Si-Ge  alloys,  A,  B  and  D  are  positive  in  sign,  but  C 
is  negative  in  sign.  Of  course,  the  values  of  A,  B,  C  and  D  change 
somewhat  with  the  chemical  composition  of  Si-Ge  alloys  [13,14] 
or  its  fabrication  method.  The  present  D  values  obtained  from 
the  temperature  dependence  of  Z  agree  with  the  resultant  values 
of  (2 A-B-C)  (see  Table  1)  within  4%.  It  indicates  that  the  present 
approximation  is  valid  for  Si-Ge  alloys. 

The  Thomson  coefficients  p  for  the  p-  and  n-type  Si-Ge  alloys 
were  calculated  using  Eqs.  (9)  and  (10).  As  a  result,  p(Th)  at 
Th  =  1000  K  was  175  pV/K  for  the  p-type  alloy  and  -237  pV/K  for 
the  n-type  one.  In  contrast,  p{Tc)  at  Tc  =  400  K  was  70  pV/K  for 
the  p-type  alloy  and  -95  pV/K  for  the  n-type  one.  The  values  of 
p(Th)  were  much  higher  in  absolute  value  than  the  p=  129  pV/K 
estimated  theoretically  by  Ioffe  [7]  using  a  non-degenerate  model, 
while  those  of  p{Tc)  were  much  lower  in  absolute  value  than 
P  =  129  pV/K.  The  band  structure  of  Si-Ge  alloys  might,  thus, 
change  from  a  degenerate  state  to  a  non-degenerate  one  with  an 
increase  of  temperature. 


Fig.  2.  Temperature  dependences  of  (a)  |a|,  (b)  p  and  (c)  k  measured  by  us  for  the  p- 
and  n-type  Si-Ge  alloys  (Ref.  [9]). 

3.2.  Dependence  of  p/rj0  on  AT,  B  and  C 

The  dependence  of  p/p0  on  AT  was  calculated  using  Eq.  (21 )  as  a 
function  of  AT  using  B  =  1.15  x  lO  ^K1  and  0-1.87  x  lO^IC1 
for  the  p-type  Si-Ge  alloy  and  B  =  1.39  x  10-3 1<-1  and 
C=  -2.56  x  10_4I<_1  for  the  n-type  [9]  (see  Table  1),  as  shown  in 
Fig.  4.  Consequently,  p/p0  of  the  p-  and  n-type  alloys  increase  par- 
abolically  with  an  increase  of  AT  and  increased  to  1.005  for  the  p- 
type  alloy  and  to  about  1.009  at  AT  =  600  K.  As  shown  in  the  figure, 
the  degree  of  contribution  from  the  term  of  B  to  p/p0  is  about  4.6 
times  larger  than  that  from  the  term  of  C.  The  degree  of  increase 
in  p/p0  due  to  the  contribution  of  B  and  C  is  a  little  lower  than 
1%  at  AT  =  600  K  for  the  present  Si-Ge  alloys.  If  the  values  of  their 
parameters  remain  unvaried  with  AT,  however,  p/p0  would  be  en¬ 
hanced  significantly  with  further  increase  in  AT. 

Next,  the  dependence  of  p/p0  on  B  and  C  was  calculated  as  a 
function  of  X(=B  =  |C|)  at  the  mean  temperature  of  T=  700  K  using 
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T  (K) 


Fig.  3.  Temperature  dependences  of  Z  measured  for  the  p-  and  n-type  Si-Ge  alloys 
(Ref.  [9]). 


Table  1 

The  values  of  A,  B,  C  and  D  at  T  =  700  K  estimated  from  the  temperature  dependences 
of  the  thermoelectric  properties  measured  for  the  p-  and  n-type  Si-Ge  alloys  (Ref.  [9]) 


Items  (K  1) 

p-type 

n-type 

A 

8.03  x  10“4 

7.69  x  10“4 

B 

1.15  x  10“3 

1.39  x  10“3 

C 

-1.87  x  10“4 

-2.56  x  10“4 

2A-B-C 

6.43  x  10“4 

4.04  x  10“4 

D 

6.18  x  10“4 

4.11  x  10“4 

Fig.  4.  rj/rjo  calculated  as  a  function  of  AT  at  T  =  700  K  for  the  p-  and  n-type  Si-Ge 
alloys  (Ref.  [9])  using  their  analytical  values  of  B  and  C. 


Fig.  5.  rj/rjo  calculated  at  T  =  700  K  as  a  function  of  X  under  various  temperature 
conditions  of  (a)  AT  =  600  K  for  the  p-type  Si-Ge  alloy  and  of  (b)  AT  =  200,  400  and 
600  K  for  the  n-type  Si-Ge  alloy,  where  X  =  B  =  |C|. 

Eq.  (21 )  under  the  conditions  of  AT  =  600  K  for  the  p-type  Si-Ge  al¬ 
loy  and  of  AT  =  200, 400  and  600  K  for  the  n-type  one  and  plotted  as 
a  function  of  X  in  Fig.  5.  rjlrj0  tends  to  increase  linearly  with  an  in¬ 
crease  of  B  but  parabolically  with  an  increase  of  |C|,  so  that  the  de¬ 
gree  of  contribution  from  C  to  rj/rj 0  is  14  times  larger  than  that  from 
B  at  X  =  2  x  10-3  K1.  Using  r]lrj0  derived  approximately  from  the 
thermal  rate  equations,  one  estimates  the  degrees  of  contribution 
from  B  and  C  to  rjlrj0.  However,  the  present  results  would  be  almost 
unvaried  fundamentally  even  if  one  employed  rjlrj0  derived  exactly 
from  them.  With  our  Si-Ge  alloys,  the  values  of  B  are  relatively  close 
to  2  x  10-3K-1,  while  those  of  |C|  are  only  about  a  tenth  of 
2  x  10-3  K_1.  Even  at  C  =  2  x  10-3  K_1,  the  present  approximation 
is  valid  because  (CAT)2/12  in  Eq.  (6)  is  0.12  at  AT  =  600  K,  which 
is  fully  small  compared  with  1.  Indeed,  there  is  a  TE  material  having 
a  |C|  value  close  to  2  x  10-3  K_1  among  the  TE  materials  for  energy 
conversion;  for  example,  |C|  =  1.8  x  10-3  I<-1  for  CoSb3  [15].  There¬ 
fore,  the  value  assumed  here  for  C  is  not  imaginary  but  realistic. 
Interestingly,  as  is  evident  from  Eq.  (21 ),  rjlrj0  of  a  TE  element  is  en¬ 
hanced  significantly  even  when  D  =  0  K_1  (for  example,  A  =  0  K-1 
and  B(=-C  =  2  x  10-3  I<-1).  Moreover,  when  B  and  |C|  have  a  large 
value  of  2  x  10-3  K1,  rjlrj0  is  found  to  be  enhanced  by  about  10% 
at  AT  =  600  K,  as  shown  in  Fig.  5.  In  any  case,  the  sign  and  magni¬ 
tude  of  B  and  the  magnitude  of  C  are  found  to  be  very  important 
to  enhance  rj  as  high  as  possible. 

3.3.  ZeffT  obtained  by  ascribing  the  increase  in  rj/rj0  due  to  B  and  C  to 
the  increase  in  the  intrinsic  ZT 

The  effective  ZeffT  for  a  TE  element  was  calculated  by  ascribing 
the  increase  in  rj/rjo  due  to  B  and  C  to  the  increase  in  the  intrinsic 
ZT.  For  simplicity,  Eq.  (21)  is  rewritten  as 
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3-  =  \  +  y  (26) 

% 

where 


(ATfy  [B(l-y) 
4(yTc  +  Th)  2 


C2r(l+y)' 

3 


(27) 


where  y  =  (1  +  ZT)  1/2.  If  we  attribute  the  increase  in  rjlrj0  due  to  B 
and  C  to  the  increase  in  the  intrinsic  ZT,  rj  is  expressed  as 


_(l-yeff)AT 

(7effT c  +  ^h) 


(28) 


using  the  effective  yeff=(l  +Zeff T)  1/2.  Substituting  Eqs.  (22)  and 
(28)  into  Eq.  (26),  the  following  equation  is  obtained  as 


7eff  =  7  - 


(1  -y)(yrc  +  rh)Y 

2  T 


(29) 


Substituting  Eq.  (27)  into  Eq.  (29),  moreover,  the  effective  ZeffT  is  ex¬ 
pressed  by  a  similar  approximation  as 


ZeffT  =  ZT 


(AT)2  / 5(1  -  y)  C2r\ 
4 T  (2(l+y)+  3  J 


(30) 


using  ZT,  B  and  C.  This  expression  shows  that  ZeffT  is  always  larger 
than  the  intrinsic  ZT  if  B  >  0.  It  indicates  that  the  presence  of  B 
and  C  results  in  a  substantial  increase  in  ZT  if  B  >  0.  In  order  to 
investigate  the  effect  of  the  increase  in  rjlrj0  due  to  B  and  C  on  ZeffT, 
ZeffT  and  ZT/ZeffT  were  calculated  as  a  function  of  ZT  at  T=  700  K 
using  Eq.  (30)  under  the  temperature  differences  of  AT  =  200,  400 


Fig.  6.  (a)  The  effective  (ZeffT)  calculated  at  T  =  700  K  as  a  function  of  ZT  using  Eq. 
(30)  by  ascribing  the  increase  in  rj  due  to  B  =  |C|  =  2  x  10-3  /C1  to  the  increase  in  the 
intrinsic  ZT  and  (b)  the  ratio  (ZeffT)/  (ZT)  calculated  as  a  function  of  ZT  using  Eq.  (30). 


and  600  K  for  X  =  8=  |C|  =2  x  10“3  IC1.  The  effective  ZeffT  is  well  de¬ 
fined  even  at  AT  =  600  K,  because  the  present  approximation  holds 
even  for  AT  =  600  K  and  X  =  2  x  10“3  K“\  as  mentioned  earlier.  As 
shown  in  Fig.  6a,  ZeffT  increases  almost  linearly  with  an  increase 
of  ZT  and  the  gradient  of  ZeffT  for  ZT  tends  to  be  enhanced  more 
strongly  at  larger  AT.  The  ratio  of  (ZeffT)/(ZT)  also  increases  more 
significantly  with  an  increase  of  the  intrinsic  ZT  at  larger  AT,  as 
shown  in  Fig.  6b.  It  suggests  that  the  effective  ZeffT  of  a  TE  element 
is  enhanced  significantly  even  when  D  =  0  K-1  (for  example, 
A  =  0  K1  and  B  =  -C=  2  x  10-3  K-1).  It  was,  thus,  found  that  the 
sign  and  magnitude  of  the  temperature  dependences  of  TE  proper¬ 
ties  have  a  significant  influence  on  the  energy  conversion  efficiency 
t]  of  a  TE  element,  particularly  at  large  AT. 

4.  Conclusions 

The  new  thermal  rate  equations  were  constructed  by  introduc¬ 
ing  newly  the  temperature  dependences  of  p  and  k  of  TE  materials 
on  the  assumption  that  they  vary  linearly  with  temperature.  The 
term  of  linear  temperature  dependence  of  a  disappears  but  those 
of  p  and  k  remain  in  the  thermal  rate  equations.  The  energy  con¬ 
version  efficiency  y\  for  a  single  TE  element  was  first  formulated 
by  approximate  analysis,  so  that  the  formula  derived  here  for  y\ 
can  be  utilized  for  estimation  of  a  TE  element  in  which  each  TE 
property  has  a  linear  temperature  dependence,  rj/rjo  was  derived 
and  calculated  as  functions  of  AT,  ZT,  B  and  C  for  Si-Ge  alloys, 
where  B  =  (dp/dT)/p  and  C  =  (d7c/dT)/7c,  so  that  rjlrj0  was  found  to 
be  enhanced  significantly  by  the  value  of  B  rather  than  C  with  an 
increase  of  AT.  When  C  becomes  equal  to  B  in  absolute  magnitude, 
however,  the  major  part  of  the  increase  in  rjlrj0  results  from  |C|  and 
rj/rjo  is  increased  by  approximately  10%  at  AT  =  600  K.  When  we  as¬ 
cribed  this  increase  in  rjlrj0  due  to  B  and  C  to  the  increase  in  ZT,  it  is 
found  that  the  degree  of  increase  in  ZT  at  AT  =  600  K  corresponds 
to  approximately  10%  at  X  =  B=  |C|  =  2  x  10_3I<_1.  The  approxi¬ 
mate  expression  obtained  here  for  rj  reflects  sufficiently  the  effect 
of  the  temperature  dependence  of  p  and  k  on  rj.  It  was  clarified  by 
taking  the  temperature  dependences  of  p  and  k  into  the  thermal 
rate  equations  so  that  they  also  contribute  to  an  enhancement  of 
the  energy  conversion  efficiency  rj. 

Appendix  A 

Let  us  assume  that  the  Seebeck  coefficient  a,  the  electric  resis¬ 
tance  R  and  the  thermal  conduction  I<  of  a  TE  parallelepiped  in¬ 
crease  or  decrease  linearly  with  an  increase  of  T.  Here,  we  will 
show  that  (fl(Th))  or  (ft(Tc))  averaged  over  a  half  (hot  or  cold)  side 
of  a  TE  element  and  the  resultant  (K(T))  averaged  all  over  an  ele¬ 
ment  are  expressed  as  a  function  of  AT.  When  we  assume  that 
the  temperature  in  a  TE  parallelepiped  changes  linearly  with  z 
along  its  length  direction,  as  shown  in  Fig.  1,  the  temperature  Tz 
at  a  position  z  is  expressed  as 

Tz  =  T+(z-L/2)AT/L.  (A.l) 

When  the  electrical  resistivity  p  changes  linearly  with  T,  p(Tz)  at  Tz 
is  then  expressed  as 

P(TZ)  =  p(T)[l  +  (Tz  -  T)B], 

where 

1  dp(Tz)  _  1  d R(TZ) 
p(T)  dTz  R(T)  d  Tz 

Therefore,  the  average  (R(Th))  of  the  upper  (hot)  side  (corresponding 
to  half  of  a  TE  parallelepiped)  is  given  by 


(A.2) 

(A-3) 
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Using  the  relation  dlz/dz  =  AT/L,  (R{Th))  is  obtained  as 


References 


(R(Th))  =  (R(T)/2)(l+BA7/4)  (A.5) 

In  the  same  way,  the  average  R(TC)  is  also  derived  as 

(R(Tc)>  =  (R(T)/2)(1-BAT/4)  (A.6) 

(R(Th))  is  higher  by  BR(T)AT/4  than  <R(rc)>  for  B  >  0. 

Next,  if  the  thermal  conductivity  k  changes  linearly  with  tem¬ 
perature,  K(rz)  at  Tz  is  then  expressed  as 

k(Tz)  =  k(T)[1+(Tz-T)C],  (A.7) 


where 

1  dK(7z)_  1  d  K(TZ) 

k(T )  dTz  K(T)  d  Tz 


(A.8) 


When  there  is  a  temperature  difference  AT  along  a  TE  element, 
therefore,  the  resultant  ( K{T ))  of  a  TE  element  is  given  by 
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(/C(T)>  "  S  io  k(Tz)  02  ~  S  JTc  k(T2)  d Tz  °  z 
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